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We present a method to detect properties of quantum noisy channels, assuming that some a priori 
information about the form of the channel is available. The method is based on a correspondence 
with entanglement detection methods for multipartite density matrices based on witness opera- 
tors. We illustrate the method in the case of entanglement breaking channels and non separable 
random unitary channels, and show how it can be implemented experimentally by means of local 
measurements. 
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The possibility of determining properties of quantum 
communication channels or quantum devices is of great 
importance in order to be able to design and operate the 
channel at the best of its performances. In many real- 
istic implementations some a priori information on the 
form of a quantum channel, or a quantum noise process, 
is available and it is of large interest to determine ex- 
perimentally whether or not the channel has a certain 
property. The aim of this work is to propose efficient 
methods to detect this possibility by avoiding full quan- 
tum process tomography, which allows a complete recon- 
struction of the channel but it requires a large number of 
measurement settings. At the same time, from the point 
of view of implementations, our procedure is experimen- 
tally feasible with present day technology based on local 
measurements . 

Quantum channels, and in general quantum noise pro- 
cesses, are described by completely positive and trace 
preserving (CPT) maps which can be expressed in 
the Kraus form IT] as 



(i) 



where p is the density operator of the quantum system 
on which the channel acts and the Kraus operators {A^} 
fulfil the constraint ^2 k A\Ak = 1. 

In order to develop the detection method proposed, we 
will use the Choi-Jamolkowski isomorphism, which gives 
a one-to-one correspondence between CPT maps acting 
on T>{K) (the set of density operators on H) and bipartite 
density operators Cj( on "H ® H. This isomorphism can 
be described as 



(2) 



where ^ is the identity map, and \a) is the maximally en- 
tangled states with respect to the bipartite space T-L (g> H, 
i.e. \a) = X}fe=i 1^)1^) ( we consider here quantum 
channels acting on systems with finite dimension d) . This 
is schematically depicted in Fig. [T] 

In the Kraus operators description the above relation 



l«> 













Figure 1. Scheme showing how the Choi-Jamolkowski iso- 
morphism works, on the left the map on the right the 
corresponding Choi state 



takes the form 



5^(^®l)|a)(a|(4®l) 



(3) 



In this work, by the above isomorphism, we link some 
specific properties of quantum channels to properties of 
the corresponding Choi states Cj(. Therefore we start 
from state detection methods to construct methods for 
quantum channels detection. In particular, we find a con- 
nection between quantum channel properties and (mul- 
tipartite) entanglement properties of the corresponding 
Choi states. We will consider properties that are based 
on a convex structure of the quantum channels. 

Consider as a first case the class of entanglement break- 
ing channels [2]. A possible definition for an entangle- 
ment breaking channel is based on the separability of its 
Choi state. More explicitly, a quantum channel is entan- 
glement breaking if and only if its Choi state is separa- 
ble. The set of Choi states corresponding to entangle- 
ment breaking channels therefore contains only bipartite 
separable states. This allows to formulate a method to 
detect whether a quantum channel is not entanglement 
breaking by exploiting entanglement detection methods 
designed for bipartite systems [3]. To this end, we re- 
mind the concept of entanglement detection via witness 
operators @]: a state p is entangled if and only if there 
exists a hermitian operator W such that Tr [Wp] < and 
TrfM^psep] > for all separable states. As a simple ex- 
ample of quantum channel detection consider the case of 
two-dimensional systems and the single qubit depolaris- 
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ing channel, denned as 



[p} = 



i=0 



(4) 



where oo is the identity operator, {<Ji\ (i = 1, 2, 3) are the 
three Pauli operators a x ,<r y , a z respectively (for brevity 
of notation in the following the Pauli operators will be 
denoted by X, Y and Z), and pa = 1 — p (with p € 
[0,1]), while Pi = p/3 for i = 1,2,3. Such a channel is 
entanglement breaking for p > 1/2. The corresponding 
set of Choi bipartite density operators is given by the 
Werner states 

P P = {l-\p)\a){a\ + P -l . (5) 

It is then possible to detect whether a depolarising 
channel is not entanglement breaking by exploiting an 
entanglement witness operator for the above set of states 
[31 |S] , which has the form 

Web = -(1®1+X®X-Y®Y-Z®Z) . (6) 

The method can then be implemented by preparing a 
two-qubit state in the maximally entangled state |a), 
then operating with the quantum channel to be detected 
on one of the two qubits and measuring the operator 
Web at the end. As we can see from the above form, the 
operator (|6| can be measured by performing a set of three 
local measurements performed on the two qubits. If the 
resulting average value is negative, we can then conclude 
that the channel under consideration is not entanglement 
breaking. 

We will now consider a different scenario, namely the 
case of random unitary channels (RU) . These are defined 



k 



(7) 



where Uk are unitary operators and pk > with ^2 k Pk = 
1. Notice that this kind of maps includes several inter- 
esting models of quantum noisy channels, such as the al- 
ready mentioned depolarising channel or the phase damp- 
ing channel and the bit flip channel [6] . Random unitaries 
were also studied extensively and characterised in Ref. 

m- 

Let us now assume that the system on which the ran- 
dom unitary channel acts is a bipartite system pab (com- 
posed of systems A and B) . We can then identify a class 
of random unitary maps which is separable, namely that 
can be written in the form 

y[pA B ] = ^2pk(Vk,A ® W KB )p(Vl A ® Wis), (8) 

k 
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Figure 2. Scheme showing how the Choi-Jamolkowski isomor- 
phism works in the case of four-partite states. The state \a) 
on the right is the maximally entangled state with respect to 
the bipartition AB — CD. 



where pab is a bipartite system, and both Vk,A and Wk,B 
are unitary operators for all fc's, acting on systems A and 
B respectively. Quantum channels of the above form are 
named separable random unitaries (SRU) and they form 
a convex subset in the set of all CPT maps acting on 
bipartite systems. Interesting examples of channels of 
this form are given by Pauli memory channels [8]. 

When considering quantum channels acting on bipar- 
tite systems, the Choi state is a four-partite state (com- 
posed of systems A, B, C and D), as shown in Fig. [2] 
Notice that the state |a) = ^= J2kj=i %3)ab% ])cd 
(where cLab = A-aAb is now the dimension of the Hilbert 



space of the bipartite system AB), can also be written as 
| a) = | a) ac I o) bd, namely it is a biseparable state for 
the partition AC — BD of the global four-partite system. 
The Choi states corresponding to SRU channels there- 
fore form a convex set, which is a subset of all bisepara- 
ble states for the partition AC — BD. Since the extremal 
maps of the set of separable random unitaries are given 
by local unitaries Ua®Ub, the extremal bipartite pure 
states in the corresponding set of Choi states have the 
form (we name this set of four-partite density operators 
Ssru) 



(U A ® 1c)\oc)ac ® (U B ® 1d)\oc)bd 



(9) 



We can now detect non separable RU maps (which cor- 
respond to Choi states that are entangled in the bipar- 
tition AC — BD) by designing suitable witness operators 
that detect the corresponding Choi state with respect to 
biseparable states belonging to Ssru i n t ne bipartition 
AC— BD. 

We illustrate this procedure with a simple example. 
Consider the case of detecting a CNOT operation, which 
is a non separable unitary operation acting on two qubits. 
The corresponding Choi state has the form 

Ccnot = |CNOT)(CNOT| (10) 
= (CNOT (g) l)|a)(a|(CNOT ® 1), 

where the CNOT operation is given by 

A 0~ 



CNOT = 



X 



(11) 
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with 1 representing the 2x2 identity matrix, and X the 
usual Pauli operator. A suitable detection operator for 
the CNOT gate can be constructed as 



CNOT 



= pi -a 



CNOT 



(12) 



where the coefficient f3 is the squared overlap between 
the closest biseparable state in the set Ssru and the 
entangled state Ccnot, namely 



P = max (</>|C C not| 

WeSsRu 



(13) 



Notice that, since the maximum of a linear function over 
a convex set is always achieved on the extremal points, 
the maximum above can be always calculated by max- 
imising over the pure biseparable states (J9J. Further- 
more, the coefficient f3 used above is the smallest value 
of the coefficient leading to a semi-positive expectation 
value of Wcnot over all biseparable states, so VFcnot is 
clearly optimal among the witnesses of this form. 

The task is now to compute the coefficient ft. This task 
can be accomplished as follows, by exploiting techniques 
to calculate the amount of entanglement in bipartite sys- 



tems. The state ( 10 ) is not separable with respect to the 
split AC — BD and it can be expressed in the Schmidt 
decomposition regarding that split as 



|CNOT) = —(\00) AC \ a ) CD + |11) ac\4> + )cd) 



(14) 



where |i/> + ) = ^(|01) + 110)). The above expression nat- 
urally proves that the maximum overlap with any bisep- 
arable state w.r.t. AC — BD cannot exceed the value of 
1/2. Since the convex set Ssru of allowed states in our 
optimisation problem is smaller than the set of all bisep- 
arable states, this would give us only an upper bound 
for the maximum overlap /?. However, two local unitary 
operations U a and Ub that saturate this bound can be 
explicitly found, namely 



U A = S, 



U b 



(15) 
(16) 



where S is the phase gate given by S = diag(l, i). This fi- 
nally proves that the optimal coefficient j3 equals 1 /2 also 
if we restrict to the set of biseparable states Ssru- There- 
fore, the detection operator Wqnot can be decomposed 
into a linear combination of local operators as follows 



Wcnot = 7rr(31 1 1 1 1 - 1 X IX - XXX 1 —X 1 XX 

64 

- ZZ1Z + ZY1Y + YYXZ + YZXY 
-Z1Z1 -ZXZX + YXY 1 +Y 1 YX 

- 1 ZZZ + 1 YZY + XYYZ + XZYY) , 

(17) 



\a) 



1 


i 

B 






5 

c 




D 







Figure 3. Experimental scheme implementing the detection 
of the CNOT gate. 



where for simplicity of notation the tensor product 
symbol has been omitted. As we can see from 
the above form, the CNOT can be detected by us- 
ing nine different local measurements settings, namely 
{XX XX, ZZZZ, ZYZY, YXYX, YYXZ, YZXY, 
ZXZX, XYYZ, XZYY}. Following it can be also 

easily proved that the above form is optimal in the sense 
that it involves the smallest number of measurement set- 
tings. From an experimental point of view, the opti- 
mal detection procedure then works as follows: prepare 
a four-partite qubit system in state \a) = \a) Ac\ a ) bd, 
input qubits A and B to the quantum channel and fi- 
nally perform the set of nine local measurements reported 



above in order to measure the operator ( 17 ). If the result- 
ing average value is negative then the quantum channel 
is detected as a non separable random unitary map. The 
experimental scheme is shown in Fig. [3j 

Notice that the number of measurements needed in 
this scheme is much smaller than the one required for 
complete quantum process tomography [6], that scales 
as d\ B . The number of measurement settings in the de- 
tection scheme can be further decreased if we allow a non 
optimal detection operator, in the sense that the coeffi- 
cient P in ( 12 ) is smaller than the maximum value. In 



this case, since the state Ccnot is a stabilizer state with 
generators {XXX 1, 1 X 1 X, Z 1 Z 1, ZZ 1 Z}, an alter- 
native detection operator can be derived, following the 
approach of Ref. [TO]. The resulting suboptimal detec- 
tion operator turns out to be 



Wcnot = 31-2 



(1+XXX1) (1 + 1 A 1 A) 



2 2 
(1+Z1Z1) (1+ZZ1Z) 



(18) 



which requires only the two local measurement settings 
{AAAA, ZZZZ}. 

We will now study the robustness of the method in the 
presence of additional noise, which can influence the op- 
eration of the quantum channel. Let us first assume that 
phase damping noise is present, acting independently on 
the two qubits A and B in general both before and af- 
ter the operation of the CNOT gate, as shown in Fig. 
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Figure 4. CNOT gate in the presence of dephasing noise. 

|1J Phase damping noise is described by a CPT map of 
the form Q, with probabilities po — q, P\ = P2 = 
and P3 = 1 — q. Notice that the global resulting channel 
shown in Fig. [4] is still a random unitary channel. 

In order to quantify the noise robustness of the wit- 
ness Wcnot with respect to phase damping noise, we 



calculate the expectation value of Wqnot given by ( 12 ) 
(with /? = 1/2) with respect to the bipartite state C^ p 
depicted in Fig. [4j i.e. the Choi state corresponding to 
the composite map M P = {9 s % <g> ^ 2 )CNOT(^ l i ® ^> 1 ). 
The problem thus reduces to evaluate the overlap be- 
tween the Choi states Ccnot and Cj( p . This can be 
easily achieved as follows. In general, given two maps ^ 
and Jzf acting on a <i-dimensional system and described 
by sets of Kraus operators {Ak} and {Bi} respectively, 
the overlap between the two corresponding Choi states 
takes the simple form 



Tr[C^C^] = i^|Tr[4B, 



(19) 



k,l 



where the double summation is over the Kraus opera- 
tors. In our case, the two maps correspond to the CNOT 
gate and the composite map defined above. After a 
lengthy calculation, this procedure leads to 



Tr[WcNOTC^ p ] = --[qjql+(l-qi){l-q2)(qi+q2-qiq2)}- 

(20) 

From the above expression we can see that 
Tr[W / cNOT'C'^] < for certain intervals of the noise 
parameters q\ and When the dephasing channels in- 
troduce the same level of noise (q\ = qi = q) the ex- 
pectation value of W-cnot turns our to be negative for 
q > 0.83 and therefore the CNOT operation can be de- 
tected in this interval. As we can see, the presence of lo- 
cal dephasing noise affects the CNOT operation in such 
a way that, above a certain amount of noise, the noisy 
CNOT operation becomes separable. 

We have also evaluated noise robustness for the am- 
plitude damping channel, which is not a random unitary 
noise and it is described by the following Kraus operators 



1 







(21) 



where 7 is the parameter characterising the damping. In 
this case, following the same procedure described above 



by considering now the composite map j&a = (3/2 
.e^CNOT^i ® ^l), we have 



Ty[W cnot C^ a ] = 



1 



16 



[(1 + \/7i72(l + V7i + V12)) 2 + 71717272] ,( 22 ) 



where we have defined 7=1 — 7. For the particular 
case of 71 = 72 = 7 we have that the above expression is 
negative for 7 < 0.31, and therefore the composite map 
can be detected as a non separable random unitary in 
this range of noise parameter 7. 

In conclusion, we have presented an experimentally 
feasible method to detect specific properties of noisy 
quantum channels, including quantum noise processes. 
The proposed procedure works when some a priori knowl- 
edge on the quantum channel is available and is based 
on a link to detection methods for entanglement prop- 
erties of multipartite quantum states via witness oper- 
ators. In this work the method has been explicitly il- 
lustrated to detect entanglement breaking or separable 
random unitary properties of quantum channels. The 
advantage over standard quantum process tomography 
is that a much smaller number of measurement settings 
is needed in an experimental implementation. Moreover, 
we want to point out that the proposed scheme relies 
on local measurements and it is achievable with current 
technology, for example in quantum optical implementa- 
tions [TT1. 
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